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The QCDSF and UKQCD collaborations have generated ensembles of gauge field configura- 
tions using Nf = 2 non-perturbatively 0{a) improved Wilson quarks and Wilson's plaquette action 
for the gauge fields. Here we want to discuss the results obtained for the axial charge of the nu- 
cleon, gA- Their interpretation is not straightforward because the quark masses in the simulations 
are larger than in nature, the volumes are somewhat smaller than infinity, the lattice spacings are 
larger than etc. So we need some guidance for the extrapolations towards the physical quark 
masses, the thermodynamic and continuum limits. Such guidance is provided by chiral effective 
field theory (ChEFT), which for selected quantities, e.g. for gA, yields parameterisations of the 
dependence on the quark mass and the volume which take into account the constraints imposed by 
(spontaneously broken) chiral symmetry. The dependence on the lattice spacing a can be included, 
but we shall not consider this possibility here. So we do not yet attempt to cope with the lattice 
artefacts remaining even after 0{a) improvement. 

If ChEFT can be successfully applied, we gain control over the chiral extrapolation and the 
approach to the thermodynamic limit. At the same time we can determine not only the physical 
value of the quantity of interest, gA in our case, but also some effective coupling constants. These 
may occur in the ChEFT expressions for other observables and be of phenomenological interest 
there. Establishing the link between Monte Carlo results and ChEFT will thus enable us to extract 
considerably more information from our simulations than just the physical value of the quantity 
under study. 

In its standard form, ChEFT describes low-energy QCD by means of an effective field theory 
based on effective pion, nucleon, . . . fields. Since the effective Lagrangian does not depend on the 
volume, besides the quark-mass dependence the very same Lagrangian governs also the volume 
dependence, and finite size effects can be calculated by evaluating the theory in a finite (spatial) 
volume. Thus the finite volume does not introduce any new parameters and the study of the finite 
size effects yields an additional handle on the coupling constants of ChEFT. The effective descrip- 
tion will break down if the box length L becomes too small, just as it fails for pion masses that are 
too large. 

The simulation parameters are listed in Table |l]. Note that we have two groups of three ensem- 
bles each which differ only in the volume. 

We compute gA from forward proton matrix elements of the flavour-nonsinglet axial vector 
current Ay^ = uy^Ysu - dy^y^d: 



The required bare matrix elements are extracted from ratios of 3-point functions over 2-point func- 
tions in the standard fashion. Compared to the computation of hadron masses, additional difficul- 
ties arise in the calculation of nucleon matrix elements such as gA '■ In general there are quark-line 
disconnected contributions, which are hard to evaluate, the operators must be improved and renor- 
malised etc. Fortunately, in the limit of exact isospin invariance, which is taken in our simulations, 
all disconnected contributions cancel in gA, because it is a flavour-nonsinglet quantity. The im- 
proved axial vector current is given by 



{p,s\Ay>\p,s\ 



(1) 



= q{x)yiiy5q{x) +acAd^iq{x)ysq{x) 



(2) 
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Table 1: Simulation parameters. 



and hence the improvement term, i.e. the term proportional to ca, does not contribute in forward 
matrix elements. The renormalised improved axial vector current can be written as 

A^=Z4(l+Z7^am)A7 (3) 

with the bare quark mass m = (1/ K"sea — 1/jQ-) / {2a). 

While the coefficient bA will be computed in tadpole improved one-loop perturbation theory, 
we calculate the renormalisation factor Za non-perturbatively by means of the Rome-Southampton 
method ^. Thus Za is first obtained in the so-called RT-MOM scheme. Using continuum 
perturbation theory we switch to the MS scheme. For sufficiently large renormalisation scales 
pL, Za should then be independent of pL. However, unless }X <^\/a lattice artefacts may spoil 
this behaviour. Since our scales do not always satisfy this criterion, we try to correct for this 
mismatch by subtracting the lattice artefacts perturbatively with the help of boosted one-loop lattice 
perturbation theory. Some lattice artefacts still remain, but we can nevertheless estimate Za- In 
Table § we compare our results with a recent determination of Za by the ALPHA collaboration [^]. 
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Table 2: Values of Za from this work and from the ALPHA collaboration. 

Our results for gA are plotted in Fig. [l|. Here nij^ has been taken from the largest available lattice 
at each (j8, K"sea) combination. The scale has been set by means of the force parameter ro with tq = 
0.467 fm, and vq/q has been taken at the given quark mass. Obviously there are considerable finite 
size effects. A similar volume dependence has already been observed in quenched simulations 

In order to describe or fit these data we use ChEFT. More specifically, we employ the so- 
called small scale expansion (SSE) [Q], which is one possibility to include explicit A(1232) degrees 
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Figure 1: Results for g^- The red triangles correspond to a spatial box length L w 1 .Ofm, the blue triangles 
belong to L « 1.3 fm, and for the black symbols the volumes are larger. The red star represents the physical 
point. 

of freedom in ChEFT. The small expansion parameter in the SSE is called e, and in 0(£^) the mass 
dependence of gA is given for infinite volume by [^] 
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where denotes the chiral limit value of gA- This expression depends on several coupling con- 
stants, all referring to the chiral limit: is the pion decay constant with the physical value of about 
92.4 MeV, Aq denotes the real part of the NA mass splitting, ca and gi are NA and AA axial coupling 
constants, respectively. Finally, C™(A) is a counterterm at the renormalisation scale A, which can 
be expressed in terms of the more conventional heavy-baryon couplings B^^{X) and Bjol'^)' 

C™(A)=B5(A)-2g«B^o(A). (5) 

Evaluating the underlying ChEFT in a finite spatial volume yields an expression for the L depen- 
dence of gA [0,1]. 

Phenomenology provides some information on the parameters appearing above. The analysis 
of (inelastic) n N scattering, in particular the process TiN nnN, suggests that choosing the 
physical pion mass as the scale A one has 

-2 



55 (A 



(-1.4±1.2)GeV" 



0. 



(6) 
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Figure 2: Results for with the finite size correction subtracted. The curve represents the first fit described 
in the text. The data points shown as open symbols have not been fitted. 



Therefore we set A = 0.14GeV in the following and identify C™(A = mf^') = ^(^(A = mf^'). In 
the real world one has Aq = 0.2711 GeV, and from an 0{£^) SSE analysis of the A width one finds 
ca = 1-5. At the physical pion mass we have gA = 1.267, while g° 1.2 [|6|]. Little is known about 
gi. In the SU(6) quark model one finds g\ = |l.2 = 2.16. For Fji one expects in the chiral 

hmitF;r~86.2MeV. 

Unfortunately, we cannot fit all parameters. So we fix Aq = 0.2711 GeV, ca = 1.5, Fj^ = 
86.2 MeV and fit g° , Bg, gi taking into account only pion masses below approximately 600 MeV. In 
contrast to Ref. [|] the physical point is not fitted. We find g^ = 1.15(12), = -0.71(18) GeV"^^ 
gi =2.6(8) with ;^^/dof = 4.23/3. Remarkably enough, these values are very well compatible with 
our phenomenological prejudices above. In Fig. ^ we plot the data with the finite size correction 
subtracted together with the fit curve. So, if the fit would be perfect the data points which differ 
only in the volume would collapse onto a single point. 

The uncertainty in g^ is large enough to cover the experimental point. To exemplify this 
circumstance we fix g° = 1.225, a value well within the range favoured by the above fit, and use 
only and gi as fit parameters. We find = -0.66(17) GeV"^ ^nd gi = 3.0(3) with x^/^^f = 
4.61/4. In Fig. || we plot the data without subtracting the finite size corrections. Using the results 
of the last fit we show curves not only for L = oo, but also for the values taken in the simulations 
for j8 = 5.29, K"sea = 0.1359. Of course, many more variations of the fit procedure are possible, but 
the overall pattern remains remarkably stable yielding g° ~ 1 .2, ~ —0.5 ... — 0.7 GeV^^, gi ~ 3 
and Fji ~ 90MeV, in accordance with phenomenological expectations. 
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Figure 3: Results for gA with curves for several values of L. 
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